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ABSTRACT 

We report on Chandra X-ray observations of possible- AGNs which have been correlated with Ultra- 
high Energy Cosmic Rays (UHECRs) observed by the Pierre Auger Collaboration. Combining our 
X-ray observations with optical observations, we conclude that one-third of the 21 Veron-Cetty Veron 
(VCV) galaxies correlating with UHECRs in the first Auger data-release are actually not AGNs. We 
review existing optical observations of the 20 VCV galaxies correlating with UHECRs in the second 
Auger data-release and determine that three of them are not AGNs and two are uncertain. Overall, 
of the 57 published UHECRs with |6| > 10°, 22 or 23 correlate with true AGNs using the Auger 
correlation parameters. We also measured the X-ray luminosity of ES0139-G12 to complete the 
determination of the bolometric luminosities of AGNs correlating with UHECRs in the first data- 
set. Apart from two candidate sources which require further observation, we determined bolometric 
luminosities for the candidate galaxies of the second dataset. We find that only two of the total of 69 
published UHECRs correlate with AGNs (IC5135 and IC4329a) which are powerful enough in their 
steady-state to accelerate protons to the observed energies of their correlated UHECRs. The GZK 
expectation is that ~ 45% of the sources of UHECRs above 60 EeV should be contained within the 
z < 0.018 volume defined by the Auger scan analysis, so an observed level of 30-50% correlation with 
weak AGNs is compatible with the suggestion that AGNs experience transient high-luminosity states 
during which they accelerate UHECRs. 



1. INTRODUCTION 

Identifying the sources of ultra-high energy cosmic rays 
(UHECRs) is one of the major outstanding goals in as- 
trophysics. Progress has proven difficult due to the rar- 
ity of high energy events and because the deflections of 
the cosmic rays as they travel through Galactic and ex- 
tragalactic magnetic fields mean that the cosmic ray ar- 
rival directions do not point back to their origins. The 
very most energetic cosmic rays (£>6x 10 19 eV) have 
attracted attention as a promising path forward. The 
energy loss due to the GZK mechanism means that CRs 
with such energies typically have traveled about 100 Mpc 
or less, significantly limiting the possible sources for such 
high energy particles. Furthermore, at these energies 
magnetic deflections of protons may be small enough to 
allow the identification of the progenitor type based on 
statistical associations. 

By 2007, the Pierre Auger Collaboration had com- 
piled enough events to begin drawing statistical corre- 
lations between UHECR events and possible sources, 
and reported a strong correlation between UHECRs and 
the ne arby galaxies listed in the iVeron-Cettv fc Veronl 
(2006) Catalog of Quasars and Active Galactic Nuclei 
(12 th Ed.) dThe Pierre Auger Collaborationl2007l[2008h . 
The scan procedure that was used steps through UHECR 
energy threshold, maximum angular separation, and 



maximum VCV galaxy redshift to find the values giving 
the lowest chance probability compared to an isotropic 
distribution, then evaluates the likelihood of such a cor- 
relation occurring by chance, by performing the same 
analysis on many isotropic datasets. In the following, the 
term "correlated" referring to a galaxy with respect to a 
UHECR simply means that the given galaxy falls within 
the angular and redshift limits Auger proposed based on 
applying the scan method to the original dataset using 
the VCV galaxy catalog. 

Of the 27 cosmic rays with energies above 57 x 10 18 eV 
that Auger detected before Aug. 31, 2007, twenty corre- 
late within 3.2° with VCV galaxies having z < 0.018 , 
about 75 Mpc (|The Pierre Auger Collaboration! 120081) . 
There are 21 VCV galaxies correlated with these 20 
UHECRs. (More than one galaxy can be correlated with 
a UHECR and vice versa.) Galaxy catalogs such as VCV 
are incomplete in the Galactic Plane, so a better compar- 
ison is obtained by restricting to |6| > 10° where the VCV 
catalog is more complete. With this restriction there 
are 22 U HECRs of which 19 UHECRs correlate with 20 
galaxies (jZaw eTaDHQOl). This is a much higher correla- 
tion than would be expected by chance from an isotropic 
source distribution, and higher t han can be expla ined by 
nearby galaxy clustering alone (jZaw et a Ll l2011|) . More 
recent Auger data continue to show a significant, albeit 
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less s trong, correlation with VCV galaxies (|Abreu et al.l 
[2Ol0h . 

An important question is whether the observed cor- 
relation with VCV galaxies implies that AGNs are the 
sources of some or all UHECRs. The VCV catalogue 
is a list of AGN candidates, so first it must be estab- 
lished whether the correlating galaxies (i.e., the VCV 
galaxies w i thin 3 .2° of a UHECR) are actually AGNs. 
iZaw et aT] (|2009T ) looked at existing observations of the 
galaxies in VCV that were correlated with the first set of 
UHECRs and found that only 14 of the 21 galaxies show 
unambiguous evidence of AGN activity. Three show no 
signs of AGN activity, while the other four have am- 
biguous optical spectra. A widely used technique for 
optical identification of AGNs, so-called BPT line ra- 
tios, com pares the relative stre ngths of diagnostic spec- 
tral lines (jBaldwin et all 119811 ). The BPT line ratios of 
the 4 ambiguous galaxies fa ll within the ranges classified 
by iKauffmann et all (|2003l ) as AGNs but they do not 
fall wi thin the line-ratio ranges adopted bv lKewlev et all 
(|2001f ) as indicative of an active nucleus. AGN activ- 
ity can often be obscured in the optical bands by dust 
obscuration, and the Kewley test excludes some known 
AGNs. Additionally, as many as half of AGNs that are 
selected based on radio or X-ray properties would not 
be identified by looking o nly at their BPT line ratios 
(Rcviglio & Hclfand 2006). Thus, determining whether 
these 4 ambiguous VCV galaxies have active nuclei re- 
quires observations outside of optical wavelengths. 

The VCV galaxies are not the only population of 
nearby galaxies foun d to co rrelate with the Auger UHE- 
CRs. iBerlind et al.l (|2010D performed an independent 
scan analysis between the Auger galaxies and nearby Lu- 
minous IR galaxies in the PSCz catalog and found an 
excess correlation. Restricting to \b\ > 10° where PSCz 
is complete, 13 galaxies of Lib, > 10 10 ' 5 £ sim correlate 
within 2.1° with one or more of the 22 Auger UHECR 
events. Some LIRGs contain an active nucleus with dust 
absorbing the AGN radiation and re-emitting it ther- 
mally, giving rise to large IR luminosities. This raises the 
question of wheth er the LIRGs f ound to correlate may 
actually be AGNs. IBerlind et all (|2010l ) showed that 6 of 
the 13 correlating IR galaxies are also in VCV and 5 of 
these do in fact host AGNfl The other 7 correlating IR 
galaxies lacked the observations needed to differentiate 
between star-formation and obscured nuclear activity. 

A key distinguishing feature of active galaxies is that 
accretion-driven radiation produces a broad spectrum ex- 
tending from the IR to the X-ray, known as the broad- 
band continuum. Normal and starburst galaxies, on 
the other hand, have broadened blackbody spectra due 
largely to stellar emission, which is peaked sharply in 
the UV/optical. This means the X-ray to optical flux 
ratios are considerably larger for AGNs than for nor- 
mal/starburst galaxies. This has been seen in detailed 
spectroscopic studi es of Chandra Deep Field sources 
(jBarger et all 12003) . If the active nucleus is obscured, 
the X-ray to optical ratio becomes even larger since dust 
absor bs UV and optical p hotons more readily than X- 
rays (jComastri et al.ll2003l ). This makes X-ray observa- 
tions, particularly when combined with observations in 

1 The one which is not a confirmed AGN falls within the VCV 
sample to be tested. 



the near-IR, optical or UV, a powerful tool f or identifying 
all classes of AGNs (jMaccacaro et al.lll988h . 

In this paper we use X-ray observations to determine 
whether the 4 ambiguous VCV galaxies, and the 7 in- 
determinate PSCz galaxies, have active nuclei. We ob- 
served ten of these possible UHECR source galaxies using 
the Chandra X-ray satellite, while for IC 5169 we used 
data from a recent XMM-Newton observation. 

Another interesting question is whether AGNs that 
correlate with UHECRs have any characteristic features 
which distinguish them from the AGNs that do not 
seem to correlate with UHECRs. A particularly rel- 
evant property is the luminosity of the AGN, as this 
helps to con strain the possible cosmi c ray acceleration 
mechanisms ([Farrar fc Gruz inov 200Sf). Previously, reli- 
able luminosities have been established for all but one 
of the correlating AGNs for UHECRs in the first data- 
release (H^etaDHoni), (jBerlind et all 12010). We ob- 
served the remaining AGN, ESO 139-G12, in order to 
obtain the first robust estimate of its bolometric lumi- 
nosity. 

2. OBSERVATIONS AND DATA REDUCTION 

Each target galaxy was observed on the Chandra 
ACIS-I detector for 5 kiloseconds. The observations took 
place between January and August of 2009. Data analy- 
sis was performed using the ciao data analysis software. 
An X-ray source was observed at all of the targets; the 
results are given in Table [T] 

For the Chandra observations, we used two different 
methods to estimate the X-ray energy flux. The direct 
method sums the energy received per unit time in the 
band of interest (2-10 keV) in a 2 arcsec window around 
the source using the eff2evt tool which takes into ac- 
count both the quantum efficiency and the effective area 
corrections of the satellite. The flux is typically domi- 
nated by higher energy and less frequent events, making 
this technique vulnerable to Poisson fluctuations when 
the count rates are small. The second technique con- 
siders only the total counts received from the compact 
nuclear source, ignoring the energies of the photons, and 
assumes that the spectrum follows a power law with a 
spectral index of 2. Then webPIMMS can estimate the 
corresponding flux, including Galactic extinction from 
webPIMMS and assuming no intrinsic extinction. For 
the galaxies which are diffuse sources, we report only the 
direct measurement of the flux within 2 arcsec of the 
galactic center, as an upper limit on the possible nuclear 
emission. 

We also analyzed the XMM-Newton observation of IC 
5179. Using the source imaging PPS processing provided 
by XMM, we took the count rate on the EPIC PN-CCD 
camera (0.5-12 keV) and used webPIMMS to estimate 
the corresponding 2-10 keV flux of the source. 

3. ANALYSIS AND RESULTS 

X-ray to optical flux ratios provide a useful way 
to classify galaxies according to their nuclear activity. 
iBarger et "all (|2002l) show that for AGNs, the ratio of 
the 2-10 keV flux to the R-band magnitude is typi- 
cally 1 > log(f x / /flj3> —1, while obscured AGNs have 

* l°g(fx/fn) = logio(f x ) + 5.5 + .4 * Rmag where f x is the 2-10 
keV flux 
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TABLE 1 
Table listing source properties 



Source Name 


Total Counts 


Hard Counts 


Flux 


Flux 


Rmag 


— r 7r ir — 

Log{f x /fR) 


£2-10 




^U.O-IU KC V j 


^Z-1U KC V ) 


erg 


^WCDr^llVllvlo ) 


— 1 

erg s 








cm" 2 s" 1 












1 7 


R 

V 


Q Q p 1 A 
O . O Hi- L^± 


1 A V 1 zL 

_L .4 Hi - 14 


1 9 z, 1 

1Z .(J 


-o.u 


/ .4 Hi oy 


Kirr 1 7^0.1 


1 1 


o 
L 


0.0 Hi -10 


1 A T? 1/1 
1.4 1 . - 14 


1 9 
1Z.4 


-0.4 


£ on 


NGC 1204 


13 


4 


1.3 E -14 


1.2 E -14 


13. 6 4 


-2.9 


5.8 E 39 


NGC 2989 


5 


2 


2.3 E -14 


4.9 E -15 


12. 5 1 ' 2 


-3.1 


8.8 E 39 


IC 4523 


6 


5 


6.5 E -14 


6.5 E -15 


13 z 


-2.5 


3.8 E 40 


ESO 270-G007 


14 


5 


2.2 E -14 




13 1 


-3 


8.4 E 39 


IC 5186 


5 





5.8 E -15" 




12.3 1 


-3.8 


3.4 E 39 


ESO 565-G006 


21 


8 


5.5 E -14 


2.4 E -14 


12. 7 1 


-2.7 


3.2 E 40 


NGC 7648 


5 


1 


3.0 E -15 




12.2 3 


-4.4 


9.7 E 38 


2MASX Jl 754-60 


3 


1 


2.7 E -15 








1.6 E 39 


IC 5179 


2685 c 






9.1 E -14 


11.4 1 


-2.98 


2.5 E 40 


ESO 139-G12 


1070 


666 


4.7 E -12 d 




12. 9 1 ' 2 


-0.2 


2.9 E 43 



Note. — Fluxes are for 2-10 keV. The first group of four rows are the correlating galaxies found in VCV. The following six are those 
from PSCz. The last row of this group, IC 5179, was observed by XMM-Newton. Only the direct flux determination method is used for 
diffuse sources. The final row, ESO 139-G12, is a known AGN whose X-ray luminosity had not previously been measured. The 2-10 keV 
luminosity, Li— 10 j is derived from the direct flux measurement; it is the total X-ray luminosity in the central region and thus an upper 
limit on X-ray luminosity of a possible AGN, a criterion for the ability to accelerate UHECRs. 

(a) No counts above 2 keV, flux given is for 0.5 - 2 keV band, (b) Johnson magnitude rather than Cousins; this does not affect the result, 
(c) 25ksec XMM-Newton observation; counts are from 0.5 - 12 keV. (d) Derived from fit, see Figure 1. 

References. — (H ILauberts fc Valentin] lfl989T ). (21 IDovle etaLl pOOBT ). (31 Ide Vaucouleurs Ik Longd ifTMl ). (41 1 Jones etHl pOOl l 



log(f x /fn) > 1. Normal and starburst galaxies have 
l°9(fx/ fa) < —2 ; making this a robust way of discrimi- 
nating AGNs, even obscured ones, from normal and star- 
burst galaxies (jDonlev et al.ll2005T ). 

The results are summarized in Figure [T] The X-ray to 
optical ratios are far below what is expected from AGNs, 
l°d(fx/ 1 Ib) > — 1) for ten of the ambiguous or indetermi- 
nate galaxies: the nine we observed with Chandra and 
IC 5179, for which we used an XMM-Newton observa- 
tion. This holds even when the larger of the two X-ray 
flux determinations is used to calculate the X-ray to op- 
tical flux ratios. The accuracies of the X-ray flux deter- 
minations are limited by shot noise, which at these low 
counts produces relative uncertainties of order one. Since 
our measure of nuclear activity depends on logio(f x ), 
this corresponds to an uncertainty of less than ~ 0.3 in 
log(f x / fii). The observed ratios are far from the AGN 
domain, so our conclusion is robust even with these large 
uncertainties. We conclude that none of these galaxies 
have active nuclei. 

The last of the indeterminate galaxies, 2MASX Jl 
754-60, lacks the optical observations needed for the 
logifx/fii) test, but its X-ray emission is diffuse and 
weak, and its near-IR properties are similar to the other 
galaxies. We therefore conclude that 2MASX Jl 754-60 
is unlikely to host an AGN, especially not one expected 
to be capable of UHECR acceleration. 

Obscured AGNs in X-ray surveys are often identified 
by the hardness ratio of the source. The hardness ratio is 
defined as H-S/H+S, where H is counts from 2-8 keV and 
S is counts from 0.5-2 keV. Typically, nearby AGN with 
column densities greater than Nh ~ 10 22 have positive 
hardness ratios, while unobscured sources have negative 
hardness ratios (jFiore et al.ll2000T ). Table [T] gives counts 
from 0.5-10 keV and 2-10 keV. No photons with energies 
above 8 keV were detected for any of the sources due to 
the small effective area of Chandra at high energies, so 



these H count values are applicable for estimating the 
hardness ratio. Although the hardness ratios have large 
uncertainties because of the small number of counts, they 
show no evidence of obscuration. 

4. AGN IDENTIFICATION FOR THE 2ND AUGER DATA 
RELEASE 

A second set of U HE events has been released by the 
Auger collaboration (|Abreu "eTalllMol l since the Chan- 
dra observations reported here were undertaken. There 
are 42 UHECRs in the new sample, of which 35 have 
\b\ > 10°. Twelve of these (all with |6| > 10°) correlate 
within 3.1° with 20 VCV galaxies with z < 0.018. We 
have reviewed existing observations of these galaxies to 
establish which host AGNs. We could not find published 
emission lines for two of them, but one (IC 4296) has 
a radio jet and therefore it has an active nucleus. The 
status of the other (G 1314-1532) cannot be determined, 
but this does not affect the correlation count because 
another, unambiguous AGN lies within the correlation 
radius of the same UHECR. 

We use the BPT optical AGN criteria compar- 
ing the 0[IH1/H/3 a nd N[II]/Ha emission line ratios 
(| Baldwin et al.l I1981D , to evaluate the nuclear activity 
in the other 18 c orrelating VCV galaxies, as was done in 
iZaw et all (|2009[ ) for the first Auger UHECRs. The re- 
sults are reported in Table [2] below. The classification is 
clear-cut for 16 of the 18 cases. For the galaxy M rk 945, 
the lower limit of IVeilleux fc Osterbrockl (|1987f) on the 
0[III]/H/3 line ratio is at the LINER-Seyfert boundary; 
we note that the actual ratio could easily be significantly 
larger. In any case, the Log(S[II] A6716 + A6731)/(Ha) 
ratio o f Mrk 945 allows it to be clas sified as a narrow-line 
AGN (|Veilleux fc Osterbrockl fl987f) . S[II] being a probe 
of the partially-ionized region produced by high-energy 
X-rays from the ionizing source which is a distinguish- 
ing feature of active nuclei. We classify the remaining 
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Fig. 1. — R magnitude vs. 2-8 keV flux for target galaxies (squares: galaxies from the VCV catalogue; circles: galaxies from PSCz; 
diamond: ESO 139-G12, the known AGN). The solid line is log(f x / fn) = —2, below which lie quiescent galaxies. The dashed line is 
l°9{fx/ fit) = — 1j which separates AGNs on the right from a transition region consisting of Starburst galaxies and weak AGNs. 



case, ESO 18-G09, as ambiguous because its emission 
line ratios satisfy the [K auffman n~et al.l ([2003D but not 
the lKewlev et all (|200lD criteria. 

Thus we conclude that 3 of the 20 VCV galaxies for 
the 2nd UHECR dataset do not have active nuclei, while 
the classification of 2 others cannot be determined with 
existing data. This "decorrelates" two of the 12 UHECRs 
but does not change the correlation status of the other 
three affected UHECRs because those remain correlated 
with other, established AGN. Altogether, as summarized 
in Table El 10 UHECRs with |6| > 10° correlate with 
confirmed AGNs, one does not, and one remains to be 
determined. 

We also provide an estimate of the bolometric lumi- 
nosity, where possible, for correlating AGNs. From this, 
Aboi _ a measure of the capability of the AGN to accel- 
erate a proton to the observed energy - is calculated. 
One of the AGNs correlating with a UHECR in the 2nd 
data-release, IC4329a, has a large enough bolometric lu- 
minosity to accelerate its UHECR if that is a proton, two 
are undetermined, and the rest do not. 



5. LUMINOSITY OF ESO 139-G12 

The AGN ESO 139-G12 is correlated with two UHE- 
CRs, making it a particularly interesting source candi- 
date. Previously only an uppe r limit on i ts bol ometric 
luminosity existed, derived by IZaw et all ((20 09) . This 
source is bright enough for us to determine its flux by 
fitting the spectrum using xspec. Assuming an absorbed 
power law, this yields a hard photon index of 0.721q 09 
(error denotes 90% confidence interval) and very low 
absorption: N H < 8.2 x 10 21 cm -2 at 90% confidence 
level (Figure [2]). From this fit we extract a 2-10 keV 
flux of 4.7±?;i x 10~ 12 erg cm" 2 s" 1 (error denotes 90% 
confidence interval). When combined with its redshift, 
z = 0.017 (« 74 Mpc assuming a standard cosmology 
with H = 70), we find the 2-10 keV luminosity of ESO 
139-G12: L 2 _ 10 « 3.1±{$ x 10 42 erg s~ 2 . 

The conversion factor relating £2—10 to the bolomet- 
ric lu minosity depends on the r elative activity of the 
AGN (|Vasudevan fc Fabianl 12007ft . In the absence of a 
black hole mass measurement for this galaxy, we follow 
IZaw et al~l ([20090 and take Lboi = 20 L2-10, as is typical 
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TABLE 2 

New Auger UHECRs and Correlated VCV Galaxies 



Year 


Day 


E CR 




Dcccr 


ICR 


ben 


VCV Galaxy 


r 


z 


VCVClass 


RcalClass 




£<bol 










(EeV) 


(J2000) 


(J2000) 


(.72000) 


(J2000) 




(cleg) 






S2 1 




erg s~ 1 






2008 


87 


82 


220.5 


-42.9 


-36.4 


+ 15.5 


NGC 5643 


2.12 


0.003 


S2 


1.2 


x 10 44 


(a) 


0.18 
















IC 4518A 


2.88 


0.016 


S2 


S2 2 


2.4 


x 10 44 


(b) 


0.36 


2008 


192 


55 


306.7 


-55.3 


-17.3 


-35.4 


IC 4995 


2.86 


0.016 


S2 


S2 3 


3.1 


x 10 42 


(c) 


0.01 


2008 


282 


61 


202.3 


-16.1 


-44.0 


45.8 


MCG -03.34.064 


1.75 


0.017 


Slh 


S2 or NLS1 4 


2.6 


x 10 43 


(d) 


0.07 
















G 1314-1532 


2.84 


0.013 


S? 


? 




— 




— 


2008 


362 


84 


209.6 


-31.3 


-40.7 


29.4 


IC 4329A 


2.19 


0.016 


SI. 2 


SI 5 


1.2 


x 10 45 


(e) 


1.70 


2009 


39 


70 


147.2 


-18.3 


-106.5 


26.6 


NGC 2989 


0.81 


0.013 


H2 


H2 6 








— 


2009 


51 


65 


203.7 


-33.1 


-46.7 


28.9 


IC 4296 


0.95 


0.013 


S3 


AGN 7 


3.2 


x 10 42 


(f) 


0.01 
















ESO 383-G18 


0.96 


0.013 


SI. 8 


S2 4 


4.6 


x 10 43 


(g) 


0.11 
















MCG -06.30.015 


1.22 


0.008 


SI. 5 


SI. 2 s 


8.0 


x 10 43 


(h) 


0.19 
















NGC 5253 


1.82 


0.001 


H2 


H2 9 










2009 


78 


61 


26.7 


-29.1 


-134.6 


-77.6 


NGC 613 


2.74 


0.005 


S? 


Composite 10 




*** 






2009 


80 


66 


170.1 


-27.1 


-80.9 


31.5 


IRAS 11215-2806 


1.52 


0.014 


S2 


S2 4 


4.0 


x 10 43 


(g) 


0.09 


2009 


212 


57 


122.6 


-78.5 


-68.8 


-22.8 


ESO 18-G09 


1.00 


0.017 


S2 


H2?ll 










2009 


219 


57 


29.4 


-8.6 


166.1 


-65.8 


IC 184 


1.85 


0.018 


S2 


S2 12 


4.8 


x 10 43 


(0 


0.15 
















NGC 788 


1.98 


0.013 


Slh 


S2 or LINER 13 


1.8 


x 10 44 


(j) 


0.55 
















IRAS 01475-0740 


2.21 


0.017 


Slh 


Slh 11 - 14 


2.9 x 10 43 


0.09 


2009 


304 


61 


177.7 


-5.0 


-83.8 


54.7 


MCG -01.30.041 


0.50 


0.018 


SI. 8 


S1.8 15 




*** 




*** 


2009 


326 


57 


5.4 


-5.6 


103.3 


-67.3 


IRAS 00160-0719 


1.63 


0.018 


S2 


H2 16 
























MRK 945 


2.43 


0.015 


S 


S2 17 


5.4 


x 10 43 


(k) 


0.17 



NOTE. — This tabic lists each Auger UHECR from the update paper (with E > 55 EeV) whose arrival directions arc within 3.1° of a nearby (z 
< 0.018) VCV galaxy, giving the year and thejulian day each was recorded, their energies and their positions (equatorial and galactic) in degrees 
( Abrcu ct al. 2010). Next are the correlating VCV galaxies for each, their separation from the UHECR in degrees, and their rcdshift, and VCV 
classification. The correct optical classification is then listed (taken from the literature where available); "?" indicates the optical classification 
is ambiguous. The last two columns show the bolometric luminosity for the AGN and Aboi, defined as L^ol X 10 45 X E2o — 2 , the ratio which 
compares the bolometric luminosity to the theoretical luminosity necessary for producing the associated UHECR. The bolometric luminosities were 
derived form [OIII] or 2-10 keV luminosities found in literature. 

REFERENCES. — (1) IWu et al.l d2011bf) , (2) IZaw et~aH 120091), (3) IWoo U Urrvl 120021), (4 ) IGu et al.l 1200611, (5) IShu et~aT] 112010'), (6) 
IGoncalves et al.l JT999I). (71 "Cavagnolo et alj 1201011, (8) IBetmert et al.l 120061), (9) IHo et all 119971), (101 IVeron et alj 119971) , (11) Ide Griip et"aT] 
(1992), (12) Kollatschny & Frickc (1987), (13) Winter et al. (2010), (14) Tran (2001), (15) Kollatschnv et al. (2008), (16) Vader ct al. (1993), 
(17) [V cillcux fc Ostcrbrock [ 1987 ). (a) IZhu et"aT1 120111), (b) IZaw et all 12009f), (cl IMcKernan et alj 1201011, (d) IWu et ah| 12011al) , (e) 
IBrightmao fe Nandraj ihoilll, (f) IPellegrinii hoia V (g) IGu et alj 1200611 . (h) [Vasudcvan c t~aTl 120091 1, (i) IKollatschnv fc Frickel 119871 1, (j) 
IKraemer eTaT] 1201111 , (k) Mulchacy et al. ( 19*9611 . 



for an AGN with Eddington ratio less than 0.1. We there- 
fore estimate the bolometric luminosity of ESO 139-G12 
to be Lboi ~ 6.2±?;| x 10 43 erg s" 1 . 

For completeness, Table [T] gives the 2-10 keV luminos- 
ity for all of the galaxies which were observed. With an 
appropriate conversion factor, this can be used to get an 
upper limit on possible weak nuclear activity which may 
be of interest in other studies of UHECR sources. 

6. DISCUSSION 

These Chandra observations complete the task of 
determining which galaxies found to correlate with 
UHECRs bylThe Pierre Auger Collaboration! (pOOl and 
IBerlind et alH 20101 ) have active nuclei, and of determin- 
ing the bolometric luminosities of the AGNs correlating 
with UHECRs in the first data-release. The X-ray fluxes 
of the ten unclassified correlating galaxies studied here 
all fall within the range typical of normal and starburst 
galaxies of the same optical magnitude, hence we find no 
evidence of activity in any of them. Of course, it is im- 
possible to rule out the possibility of very low luminosity 
nuclear activity which is energetically dominated by the 
host galaxy, but such weakly active AGNs are not well- 
motivated candidates for UHECR acceleration anyway 
(|Farrar k. Gruzinovl[2009l ) . 

The result, then, is that only 14 out of the original 
27 Auger UHECRs (13 out of 22 UHECRs with |6| > 
10°) correlate with an actual AGN, using the Auger scan 
parameters to correlate UHECRs with VCV galaxies but 
discarding candidate sources which are not in fact AGNs. 



Of the 13 highly luminous IR g alaxies (Lt R > lO n - 5 L ) 
found bv IBerlind et al.l (|2010f ) to correlate within 2.1° 
with one (or more) of the 22 UHECRs with |6| > 10°, 
we find that only five are also AGNs0 Thus five of the 
27 UHECRs, and one of the 22 UHECRs with |b| > 10°, 
have neither an AGN nor a LIRG within 3.2° or 2.1° 

respectively. 

From the second Auger data set (jAbreu "eTani20Tol) . 
we find 10 UHECRs with |6| > 10° correlate with con- 
firmed AGNs, one does not, and one remains to be de- 
termined. 

Thus using the full dataset of UHECRs with |6| > 10°, 
where Galactic extinction does not hide source candi- 
dates, (13 + 10)/(22 + 35) = 0.40g;|| of the UHE- 
CRs correlate with a confirmed AGN using the correla- 
tion parameters proposed in the origi nal Auger analysis 
(The Pierre Auger Collaborationll2007l ) , where the upper 
and lower valu es are 1 — a limits for Poisson statistics 
(|Gehrelsl [T986). It is interesting that when actual AGNs 
rather than VCV galaxies are used, the agreement be- 
tween the correlation found in the early and later Auger 
datasets becomes closer. Using the correct galaxy attri- 
butions found here, the first and second data sets indi- 
vidually give correlations to confirmed AGNs of 13/22 = 

°- 59 8:43 and 10 / 35 = °- 29 a20 which differ at the 1 - 07 " cr 
level, compared to the individual datasets being 1.63-(7 
from the mean using the uncorrected VCV-attribution. 

2 We could not run this correlation for the full 27 UHECRs 
because IRAS does not observe within the Galactic plane. 
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Fig. 2.— Data and fit of the spectrum of ESO 139-G12 to a power law at the source (r = 90% confidence) with N H < 



8.2 X 10 21 cm 2 at 90% confidence level. This is our best fit model of the data, with a reduced \ 2 of 1.7. This fit corresponds to a 2-10 
erg cm -2 s — 1 (error denotes 90% confidence interval). 



keV flux of 4.7±°;f x 10" 



If the last VCV galaxy is confirmed as an AGN, the cor- 
relation would rise to (13 + ll)/(22 + 35) = 0.42g;g of the 
UHECRs correlating with an AGN, and 11/35 = 0.3lg;^ 
in the second dataset alone. 

Given the uncertainty in the fraction of VCV galaxies 
which are actually AGNs, and the fact that the scan pa- 
rameters were determined prior to removing non-AGNs, 
it is difficult to assess the significance of the final cor- 
relation. Further complicating correlation studies is the 
incompleteness of source catalogs (especially within the 
Galactic plane), the composition uncertainty, and mag- 
netic deflection that can obliterate the angular correla- 
tion between the CRs and their true source, particularly 
for UHECRs with charge Z > 1. Indeed, the new Galac- 
tic magnetic field model of Jans son fc Farrarl (|2012[ ). with 
a more general form for the field constrained by exten- 
sive, all sky RM and polarized synchrotron emission data, 
predicts that Galactic deflections are small in some por- 
tions of the sky but are large in others, even for protons. 

The one established AGN which we observed in order 
to determine its bolometric luminosity, ESO 139-G12, 
proves to have Lboi compara ble to most of the other 13 
correlated AGNs examined in IZaw et all (|2009D . Knowl- 
edge of Lboi and the energies E = E 2 o 10 20 eV of the 
correlated UHECRS (89 and 59 EeV) allows us to evalu- 



ate A 



10" 



_l w 45 Lhr> ] Er,n , the figure-of- merit introduced 
by IZaw et al.l ([20091 ) to quantify the ability of an AGN 
to accelerate a proton to the energy of the correlated 
UHECR. A value of Ab ? i > 1 satisfies the accelera tion cri- 
terion for protons (c.f.. lFarrar fc Gruz inov (2009])). With 
Aboi ~ 0.1, ESO 139-G12 is thus marginal according to 
standard UHECR acceleration mechanisms for protons. 

GZK energy losses imply that (taking sources to be 
uniformly distributed in redshift) about 45% of the 
sources of protons above 60 EeV should have z < 0.018. 
Thus - taking at face value the 30-50% correlation 
we find between UHECRs and confirmed (albeit weak) 
AGNs - all UHECRs could have been produced in galax- 
ies presently hosting (generally weak) AGNs, consistent 
with the picture of transient production of UHECRs via 
exceptional, po werful flares flares in very weakly or non- 
active galaxies (jFarrar fc Gruzinovl [2009") . Examples of 
tidally produced flares in quiescent gala xies have recently 
been discovered in archival SDSS data (|van Velzen et al.l 
I2011D and obs erved (apparently in b lazar mode) by th e 
Swift satellite ([Burrows et al.ll2011[ ). (IBloom et al.l201lD 
The sp ectral energy distribution of the Ivan Velzen et al.l 
(|201lD flares are well-fit by a thin accretion disk model 
and the resultant bolometric luminosities amply satisfy 
the minimum luminosity requirement for UHECR accel- 
eration in both cases (GRF, in preparation); depend- 
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ing on how rapidly the evidence of the accretion episode 
disappears, the host galaxy of a tidal disruption flare 
may or may not show evidence of weak AGN activity in 
later observations such that it would appear in a catalog 
such as VCV. It is also possible, given the uncertain- 
ties, that other candidate sources not associated prefer- 
entially with AGNs may be responsible for some, most, 
or all UHECRs. Indeed, other studies have found corre- 
lations between the arrival directions of the Auger UHE- 
CRs and a variety of extragalactic sources: HI galaxies 
(Ghisellini et al. 2008), AGNs from the Fermi catalog 
(Nemme n et al. 2010), S wift-BAT AGNs and 2MRS 
galaxies (jAbreu et al.ll2010T ). 

7. SUMMARY 

Our observations and analysis using the Chandra 
X-ray satellite and other data establish that one- 
third of the 21 galaxies in the Veron-Cetty Veron 
catalog of AGN candidates found to correlate with 
UHECR arrival directions in the first A uger data release 
(|The Pierre Auger~ Collabor ation! [2~007[ ). do not in fact 
have active nuclei. Combining this with our measure- 
ment of the X-ray luminosity of ESO 139-G12, an AGN 
correlating with two UHECRs, implies that only one of 
the 27 UHECRs in that first Auger data-release corre- 
lates with an AGN (IC 5135) which is powerful enough 
in its steady-state to accelerate protons to the observed 
energi es, according to conventio nal acceleration mecha- 
nisms (|Farrar fc Gruzin ov 2009) . None of the correlated 



PSCz source candidates we observed have active nuclei. 

We reviewed optical observations of the correlating 
VCV galaxies for UHECRs in the second data release to 
identify which galaxies actually have an active nucleus 
- 3 do and two are ambiguous - and to determine their 
bolometric luminosities. For one of them, IC 4329a, Lboi 
is high enough to accelerate a proton to the observed en- 
ergy in the steady state; Lboi is not determined for two 
AGNs, and is insufficient in the remaining cases. 

Combining the first and second Auger data releases 
and adopting the original Auger scan parameters, our 
results are compatible with 30-50% of UHECRs being 
produced by genuine but weak AGNs within z < 0.018. 
A consistent picture is that a significant portion of UHE- 
CRs are pro duced in the flaring state of otherwise very 
weak AGNs (jFarrar fc Gruzinovl [20091 ) . since about 45% 
of the sources of UHE protons should have z < 0.018. 
However until much larger UHECR datasets and more 
complete and pure catalogs of potential source candi- 
dates become available, the source(s) of the highest en- 
ergy cosmic rays will not be settled. 
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